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Abstract. A generalisation of the odd Bernoulli polynomials related to the 
quantum Euler top is introduced and investigated. This is applied to compute 
the coefficients of the spectral polynomials for the classical Lame operator. 

1. Introduction 

The classical Bernoulli polynomials can be denned through the generating func- 
tion 

ze ** _^B k {x) zh _ 



e z - 1 ^ fe! 

fc=0 

1 1 3x 2 x 
B (x) = l, B 1 (x)=x-~, B 2 (x)=x 2 ~x+-, B 3 (x)=x 3 ^ + ^- 

(see e.g. PED- They appear naturally in the calculation of sums of powers of the 
natural numbers S k -i(n) = l fc_1 + 2 fe ~ 1 + • • • + n k ^ 1 in a simple way: 

S k - X {n) = {B k {n + l)-B k )/k, 

where B k = B k (0) are the Bernoulli numbers. All odd Bernoulli numbers except 
B\ = — ~ are known to be zero, so the odd Bernoulli polynomials (up to a multiple 
k) can be thought of as an "analytic continuation" of the sums of powers from 
natural argument n to real (or complex) x. 

In this paper we introduce a new class of polynomials, which can be considered 
as an elliptic generalisation of the odd Bernoulli polynomials B 2 k+i(x). They are 
related to the quantum top and to the classical Lame operator 

d 2 

Ls = —^+s(s + l)p{z). 
where p is the Weierstrass elliptic function 0, satisfying the differential equation 

(p') 2 = 4p 3 - g 2 p - g 3 - 

It is well-known after Ince ^U] that the Lame operator (considered on the real 
line shifted by the imaginary half-period) for integer s has a remarkable property: 
its spectrum has exactly s gaps. The ends of the spectrum Ej correspond to the 
doubly periodic solutions of the Lame equation (so-called Lame functions) . We will 
call the corresponding polynomials R,2 S +i{E) = Y[jLo(E — Ej(s)) Lame spectral 
polynomials. The computation of the polynomials J?2s+i(-E) for given s = 1, 2, 3, .... 
goes back to Hermite and Halphen j2H]- In more recent time this was investigated 
within the finite-gap theory initiated by S.P. Novikov's work ^1 (see ESI 121] for 
the latest results in this direction). 
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Here we consider a related but different problem: we would like to express the 
coefficient b k of the spectral polynomial R 2s+ i(E) = E 2s+1 + biE 2s + b 2 E 2s - 1 + 
• • ■ + &2s+i as a function of s (and thus for all values of parameter s). We will show 
that in this relation naturally appear some new polynomials, generalising the odd 
Bernoulli polynomials. 

The following remarkable relation between the Lame equation and the quantum 
Euler top, going back to Kramers and Ittmann |12| . will be crucial for us. Consider 
the quantum mechanical Hamiltonian of the Euler top (see e.g. 15 ) 

H = aiMj + a 2 Mf + a 3 M|, 

where the angular momentum operators Mj satisfy the standard commutation re- 
lations [Mi,M 2 ] = iM 3 , [M 2 ,M 3 ] = iMi,[M 3 ,Mi\ = iM 2 (we assume H = 1 for 
simplicity). 

The operator H naturally acts in any representation of the Lie algebra so(3). 
In particular, it acts in the representation space with spin s of dimension 2s + 1 
as a finite-dimensional operator H s . The claim is that if the parameters a, = e,; 
are the roots e±, e 2 , e 3 of the equation 4p 3 — g 2 p — 33 = 0, then the characteristic 
polynomial of the operator H s coincides with the spectral Lame polynomial: 

(1) det(A/ - H s ) = R 2s+1 (X). 

We discuss this in more detail in the next section. 

The Weierstrass condition e\ + e 2 + £3 = is unnatural from this point of 
view (and moreover contradicts the "physical" condition of positivity of a,), so 
we consider the case when the parameters ai are arbitrary. Let us introduce new 
parameters gi,g 2 ,g3, which are symmetric functions of ai,a 2 ,a 3 defined by the 
relation 

(2) 4(z - ai)(z - a 2 )(z - a 3 ) = 4z 3 - giz 2 - g 2 z - g 3 . 

We define the elliptic Bernoulli polynomials B 2 k+i as the coefficients in the expan- 
sion of the trace of the resolvent of H s at infinity: 

(3) tr(XI - H s yi = ± 

k=0 

or, equivalently by the relation 

B 2k+1 {s) = trH k s . 

E> 2 k+i is a polynomial in s of degree 2k + 1 with the coefficients, which are poly- 
nomials in 31,32,33 with rational coefficients. Strictly speaking we should write 
B 2 k+i(s] 31, 32, 33) rather than B 2 k+i(s), but we will use both notations depending 
on the context. When 32 = 33 = these polynomials reduce up to a factor to the 
classical odd Bernoulli polynomials: 

S 2fe +i(s;3i,0,0) = (2fc + jv^fcZl B 2k +i(s + l). 

The corresponding elliptic curve T given by the equation 

y 2 = 4x 3 - gix 2 - 3 2 ir - 33 

degenerates to a rational curve in this case. If 31 = we have the standard Weier- 
strass form of an elliptic curve. The polynomials B 2 k+i (s; 0, 32, 33) are called reduced 
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elliptic Bernoulli polynomials and denoted as £>^ +1 (s; 52, 33) (W is for Weierstrass) : 
Bf = 28 + 1, Bf =0, 

= ±-g 2 s{s + l)(2s - l)(2s + l)(2s + 3), 

K = ^53s(s + l)(2a - 3)(2a - l)(2a + l)(2a + 3)(2a + 5), 

B ™ = T^92 s ( s + !)( 2s ~ !)( 2s + !)( 2s + 3 )( 4s4 + 8s3 - lls2 - 15s + 21 )- 
1680 

The coefficients of B^+i are homogeneous polynomials in 32,33 of weight 2k if 
we assume as usual that the weights of 32 and 33 are 4 and 6 respectively (in other 
words, they are modular forms of weight 2k, see e.g. QHl)- Two interesting special 
cases 32 = and 33 = are called lemniscatic and equianharmonic respectively and 
correspond to elliptic curves with additional symmetries. 

We will present some effective ways to compute the elliptic Bernoulli polynomials, 
investigate their properties and then apply them to the calculation of the coefficients 
of the Lame spectral polynomials. In particular we prove that the coefficient bk = 
bk(s) of the Lame spectral polynomial B,2 S +i(E) = Yl^ (E - Ej(s)) = E 2s+1 + 
biE 2s + &2-E 125-1 + ••• + &2s+i is a polynomial in s, 32, 33 with rational coefficients. It 
can be computed using the reduced elliptic Bernoulli polynomials by the following 
recurrence relation with bo = 1: 

k 



The first coefficients are 61 = 0, b 2 = -^s(s + l)(2s - l)(2s + l)(2s + 3), 

, 93 

63 = 1 

d 840 



93 -s(s + l)(2s - 3)(2s - l)(2s + l)(2s + 3)(2s + 5), 



64 = 7- s(s - 1) (s + 1) (2s -1) (2s + 1) (2s + 3) (56s 4 + 76s 3 -94s 2 + 201s + 630) 

(see more of them below in section 5). 

Note that once the coefficients bk(s) are known for k = 0, 1, 2s one can find 
the eigenvalues of the quantum Euler top in the representation with spin s (integer 
or half-integer) by solving the corresponding algebraic equation R2. s +i(E) = 0. 

We conclude with the discussion of possible relations and further developments. 

2. Lame equation and quantum Euler top 

The observation that the Lame equation is closely related to the quantum top 
was done by Kramers and Ittmann at the early age of quantum mechanics |12j (see 
also |23| )■ They showed that the corresponding Schrodingcr equation is separable 
in the elliptic coordinate system and that the resulting differential equations are 
of Lame form. We are going here to re-derive this result and reformulate it in the 
modern terms. 

Consider the Hamiltonian 

H = a x Mf + a 2 Mi + a 3 M§ 
acting in the space of functions on the unit sphere 
(4) ql + q 2 2 +ql = h 
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using the standard representation of the angular momenta as the first order differ- 
ential operators 



Mi = 


-«0?2<9 3 


- 53<9 2 ) 


M 2 = 


-t(q 3 dx 


- qxd s ) 


M 3 = 


-i(qid 2 


- 92<9i )■ 



Let us introduce the elliptic (or sphero-conical) coordinates u\, u 2 on this sphere 
as the roots of the quadratic equation 

(5) + + = 0, 

d\ — U CL2 — u ay, — u 

where the parameters ax, a 2l a 3 are the same as in the top's Hamiltonian. One has 
then the following expressions for the cartesian coordinates in terms of iti,tt 2 : 

2 



[ax 


- Ul)(ai 


- ua) 


{ax 


- a 2 )(ai 


- (X3) 


(«2 


- ui)(a 2 


- U2) 


(02 


- ai)(a 2 


- 03) 


(«3 


- ui)(a 3 


- u 2 ) 


(«3 


- ai)(a 3 


- a 2 ) 



(6) ^ 

The system has an obvious quantum integral (Casimir) M 2 — "}2 M , which is 
the square of the total angular momentum operator: 

[M 2 ,Mi] = 0. 

One can check that that in the elliptic coordinate system the operators H and M 2 
have the form 

d , 1 — —, — r d . 1— — r d , 1— — r d 



(7) M 2 = [^ph Ji- (vQw ) + ^p^j^^p^jJL.)] 

ux - u 2 dux aux ou 2 au 2 

(8) 

H = — hv^^-fV 1 ^!)/-) + «iV%)/-(V^M/-)] 

ux — u 2 oux dux ou 2 au 2 

where P(u) = (u — <ii)(ti — a 2 )(u — a 3 ). Note that the operator M 2 corresponds to 
the standard Laplacian —A on the unit sphere. 

Since M 2 and H commute, one can look for joint eigenfunctions. The spectral 
problem M 2 ip = flip is well-known in the theory of spherical harmonics (see e.g. 
|18j). It is known that the spectrum has the form /1 = s(s + 1) for non-negative 
integer values of s. The dimension of the corresponding eigenspace V s is 2s + 1 and 
V s is an irreducible representation of dimension 2s + 1 of the rotation group SO3 
called representation with spin s. 

It turns out that the joint eigenvalue problem 

M 2 (f) = s(s + 1)0 



H(j) = E(f) 
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is separable in the elliptic coordinates ui,U2 (see ^J|2Z]). Namely, if we substi- 
tute (f>(ui,U2) = 0i(ui)</ ) 2(u2) into this system we find that each of the functions 
4>i (ui), 02 (^2) satisfies the same differential equation: 

(4[P(«)]*-f Jl) _ s ( s + !) u + E)^ = 0, 

du du 

which can be rewritten as 



d 2 1. 1 1 1 , d 1 s(s + l)u-E 

-T^V + 1 1 \T~V = 7 7 77 77 7'' 

du 2 ti-fli U — CL2 u — 0,3 du 4 (u — ai)(u — a2){u — 0,3) 



A remarkable fact is that this is an algebraic form of the following slightly gen- 
eralised version of the Lame differential equation 

(10) -^ + s( S + 1)p4z)^ = E^ 
where p*(z) is a solution of the differential equation 

(11) (p*) 2 = 4(p* ~ ai)(p* - a 2 )(p* - a 3 ). 

Indeed, after the change of variables it — p*{z) the equation (|10() coincides with 
(0 (see UHl)- When the sum oi + ct2 + = the equation determines the 
Weierstrass elliptic function p(z), otherwise it differs from it by adding a constant. 

It is well-known (see e.g. [5]) that for ^ to be a regular solution on the sphere 
the corresponding ip must be doubly-periodic, which implies that s is integer and 
E must have one of the 2s + 1 characteristic values E m (s). For each E m (s) there 
exists exactly one (up to a factor) doubly-periodic solution to the Lame equation 
£™(m), which is called the Lame function. Therefore the basis of the eigenfunctions 
of the operator H in the invariant subspace V s consists of 2s + 1 solutions <fr(ui, U2) 
of the form £™(ui)£™(u2)' They are called sometime ellipsoidal harmonics (see 

EE1) 

Thus, we come to the following result (cf. |121 127| ): 

Theorem 1. The characteristic polynomial of the quantum top Hamiltonian H s in 
the representation space with integer spin s coincides with the spectral polynomial 
J?2s+i(A) = Iljio(A — Ej(s)) of the generalised Lame operator \1U\) . 

Remark 1. Turbiner [25| has discovered a similar but different relation of the 
Lame equation with certain quadratic elements of the universal enveloping si (2). 
The Lame spectral polynomials are known to be factorisable and Turbiner's result 
gives an interesting interpretation for the factors in these terms. 

Remark 2. A simple relation between the quantum Euler top and the Lame 
equation mentioned above is a bit misleading. Indeed there are several spectral 
problems related to the Lame equation. We have considered only smooth real 
periodic version related to real x shifted by the imaginary half-period. If we would 
consider x just real, we would have a singular version (since p has poles on the real 
line), whose spectrum has nothing to do with the quantum top. In its turn, the 
quantum Euler top in the representation with half-integer spin s has eigenvalues 
which are just some special double eigenvalues of the periodic Lame operator, which 
in this case has infinitely many gaps. 
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3. Elliptic Bernoulli polynomials 

We define now the elliptic Bernoulli polynomials B2k+i(s) as the traces of the 
powers of H s , where H s is as before the quantum top operator H in the represen- 
tation with spin s: 

(12) B 2k+l {s; gi ,g 2 ,g z )=trH k , k = 0,1,2,... 

Here the parameters g\ = 4(01+02 + 03), g 2 = —4(0102 + 0203 + 0103), gs — 4010203 
are defined by the relation 

Theorem 2. The trace trH k is a polynomial in s of degree 2fc+ 1 anti- symmetric 
with respect to s = — ~ } whose coefficients are polynomials in <?i,<?2,33 with ratio- 
nal coefficients. When g 2 = <?3 = it reduces (up to a factor and shift) to the 
corresponding classical odd Bernoulli polynomial: 

k 

(13) B 2k+1 (s- 5 i, 0, 0) = — -A__. a afc+1 ( fl + 1). 

The first part essentially follows from the Harish-Chandra general results |5| (see 
also 0], page 268), but we give here a direct proof. 

Consider the standard basis in V s consisting of the eigenvectors \j > of M3: 
M$\j >= j\j >, j — — s, — s + 1, . . . s — 1, s. In this basis, the Hamiltonian if 
is a tri-diagonal symmetric matrix H = H s with the following elements (see e.g. 
Landau-Lifshitz ^S], page 417): 

(14) < j\H\j >= i(oi + a 2 )[s(s + 1) - j 2 ] + o 3 j 2 

< j\H\j+2 >=< j+2\H\j >= ^(01-02) V(s - j)(s - j - l)(s + j + l)( s + j + 2). 

Note that both expressions are symmetric with respect to s = — |; they are also 
homogeneous polynomials of degree 1 in 01,02,03. Now, consider any diagonal 
element of H k ; it has the form: 

<j\H k \j>= ]T < j\H\i x ><i x \H\i 2 > <i k _ 1 \H\j> 

ii,i2,---,ih-i 

where the distance between 2 consecutive indices ii, i;+i is either or ±2. Since the 
starting point and the ending point coincide, if the matrix element < ii\H\ii + 2 > 
appears along the path so does the element < i\ + 2\H\ii > . This proves that 
the diagonal matrix elements of H k are polynomials of degree 2k in both s and j. 
From 114fl they are symmetric with respect to s — — g and homogeneous symmetric 
polynomials of degree k in a±, a 2l 03. Now summing over j = — s, —s + 1, . . . s — 1, s 
and taking into account that the sums of the odd powers of j are zero while the sums 
of even powers 21 are the odd Bernoulli polynomials B 2 i+\(s + 1) (multiplied by 
2l 2 +1 ) we have the first statement of the theorem. The anti-symmetry of B 2 k+i(s) 
with respect to s — ~ | follows from the well-known property of the Bernoulli 
polynomials: B m (l — s) = (—l) m B m (s). 

In the case when ai = 02 = 0, we have g 2 — gz — 0, g\ — 4a3 and H = a^M 2 . 
The spectrum of H s is then very simple: Xj — a^j 2 for j = — s, — s + 1, . . . , s — 1, s. 

Since the sum Ylj=i i 2k ~ 2 k+i ^ 2k + 1 ^ s we ^ nus orj l am dU- This completes 
the proof of Theorem 1. 
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Note that from the point of view of the elliptic curve T given by the equation 

y 2 = 4x 3 - g x x 2 - g 2 x - g 3 , 

the last case corresponds to the limit when one of the periods goes to infinity 
("trigonometric limit"). There are two more interesting special cases: lemniscatic 
when gi = g% — and equianharmonic when g\ = g 2 = 0, corresponding to the 
elliptic curves with additional symmetries. 

It is natural also to consider the Weierstrass reduction g\ = 0; we will call the 
corresponding polynomials B\^ +1 (s; 52, 53) = B 2 k+i {s; 0, g 2 , 93) the reduced elliptic 
Bernoulli polynomials. 

Theorem 3. The elliptic Bernoulli polynomial B 2 k+i ho-s the following properties: 

(1) as a polynomial in (71,32,53 B 2 k+i is homogeneous of weight 2fc, where the 
weights of g\ , g 2 and 33 are assumed to be 2, 4 and 6 respectively, 

(2) B 2 k+i for fc > 1 is divisible by s(s + l)(2s + 1), 

(3) in the reduced case BW +1 is divisible by s{s + l)(2s — l)(2s + l)(2s + 3) for 
all fc and by s(s + l)(2s - l)(2s + l)(2s + 3) (2s - 3) (2s + 5) for odd fc, 

(4) in the lemniscatic case B 2 k+i(s; 0,^2,0) = for odd integer fc, 

(5) in the equianharmonic case B 2 k+i(s; 0, 0, 33) = if fc is not divisible by 3, 

(6) in the isotropic case a\ = a 2 = 03 = a i.e g\ = \2a,g 2 = —12a 2 , g% = 4a 3 , 
B 2k+1 (s) = a k {2 S + \)s k { S + \) k . 

The proof of the first two claims follows from the definition and the anti-symmetry 
property. To prove the third one consider the representation with spin s = 5. It is 
easy to check that H acts as the 2 by 2 scalar matrix \ {a\ +a 2 + a3)Id, which is zero 
in the reduced case. Therefore B\%+i{\) — for all fc. By anti-symmetry with re- 
spect to — I we also have £?2fc+i(~ f ) = 0- For half integer s, we know from Kramers' 
theorem (see |15j . paragraph 60) that the eigenvalues are no longer distinct but are 
double roots. For the particular case s = 3/2, these eigenvalues take the values 
±-\/[3(af +a| + a 2 )/2] (see HHJ, page 419) therefore for odd fc, <B^ +1 (3/2) = and 
again by anti-symmetry S^ +1 (— 5/2) = 0. The lemniscatic and equianharmonic 
cases follow from the first claim. In the isotropic case H s = as(s + l)Id, which 
implies the last statement. 

In the general case the elliptic Bernoulli polynomials are not zero and their 
highest coefficients are described by the following 

Theorem 4. The leading term of the elliptic Bernoulli polynomial B 2 k+i{s) = 
A s 2k+1 + Ais 2k + • • • + A 2s can be written 

(15) A s 2k+1 = 2 f Res rW - f)Z + (™ 2 +/?/)+ 7 (* 2 - f)Vfdj, 
Jo 

where a = |(ai + a 2 ), (3 = i(2a 3 — ax — a 2 ), 7 = \(a\ - a 2 ). 

Indeed, for large s and j the leading behaviour of the matrix elements of H is 
< j\H\j >= ~(ax + a 2 )[s 2 - j 2 } + a 3 j 2 = as 2 + (3j 2 , 

< j\H\j + 2 >=< j + 2\H\j >= - A {ax - a 2 )(s 2 - j 2 ) = 7 (s 2 - j 2 ). 

Therefore the leading term of the diagonal element < j\H \j > coincide with the 
constant term of the Laurent polynomial [y(s 2 — j 2 )£+ {as 2 + (3j 2 ) +7(s 2 — j 2 )£,~ 1 ] k 



s 
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in auxiliary variable £. Replacing the summation over j by the integration, which 
is fine in the leading order, we come to our formula. 

Note that from this formula the fact that the final result is a symmetric function 
of ai, a 2l a 3 (and thus is a polynomial in g±, g 2l 53 ) is not obvious at all. 

Remark. From the quasi-classical arguments we can write the highest coefficient 
Aq as the following integral over the unit sphere 

(16) A = — [ H k dn = — I (aiMf + a 2 Ml + a 3 M$) k dn, 
2tt J\m\ 2 =i 27r J|A/| 2 =i 

where dfi is the area element on the unit sphere. Thus our formula l|15|l gives 
an expression for this integral. It would be interesting to compare it with the 
calculation of this integral using elliptic coordinates. 



4. Effective way to compute the elliptic Bernoulli polynomials 

Although the definition of the elliptic Bernoulli polynomials themselves gives a 
way to compute them as traces of powers of the given matrices H s , it does not seem 
to be as effective as the following procedure based on the fact that the matrix H s 
is tri-diagonal. 

Indeed, in the standard basis \j > of the space V s the eigenvalue problem Hijj = 
Xt/j leads to the following difference equation: 

(17) C n - 2 4>n-2 + V n 1p n + C n t[) n+2 = Xlpn, 

where 

c„ = ai ~ ^ yj(s - n)(s -n - l)(s + n+ l)(s + n + 2), 

v n = i(oi + a 2 ) 2 [s(s + 1) - n 2 ] + a 3 n 2 . 

For such an equation one can use the standard procedure (see e.g. from the 
theory of solitons to find the local spectral densities, which are difference analogues 
of the famous KdV densities [2D ■ In our case it works as follows. 

Let Xn = c "^" +2 , then the equation (|T7J) becomes 

(18) c£_ 2 + (V n - X)Xn-2 + XnXn-2 = 

We look for a solution in the form Xn = X — X^fclo Xn,kX~ k . Substitution of this 
expression into the equation ltTH|l gives Xn,o — v n ,x n ,i = c^-2>Xn,2 = c^_ 2 u„_ 2 , 
and for general k > 1 the recurrence relation: 

k 

(19) Xn,k+1 — y^^Xn,iXn-2,k-i- 

2=1 

Let X = JX Xn,feA-( fc+1 ) so that X n - A(l-X) and log X „ = logA-^x 4 1 - 
Thus we have 

oo _ 

(20) log X „-logA = -^^i, 

i — 1 

where X nA = f n , ^n,2 = c 2 n _ 2 + ^, = c 7 2 l _ 2 ^ n _ 2 + v n c 2 n _ 2 + ^f, 
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On the other hand one can check that Y[ n ^\ — IT m (l — Em ^ ) where the E m (s) 
are the eigenvalues of H s . Thus 

£(^-io g A) = -£E§=E# 

n n i=l i=l 

Comparing this with the equations l|2U|) . we obtain 

s 

TrH k s = k^^k = k Zn,k- 

n n—~s 

Theorem 5. The elliptic Bernoulli polynomials £>2fc+i can be computed as 

s 

(21) B 2k+ i = k £ J »,fc, 

n=— s 

where I n , k are the local densities determined by the relations \liA \2U\) . 

This gives a very effective way to compute the elliptic Bernoulli polynomials 
since the local densities are polynomials in c\ and v n (and hence in n) and thus the 
summation over n can be done with the use of the standard Bernoulli polynomials. 
We had applied this procedure to find the first 10 elliptic Bernoulli polynomials 
using Mathematica (see 8 of them in the Appendix). 

5. APPLICATION: COEFFICIENTS OF THE LAME SPECTRAL POLYNOMIALS 

We will consider again the generalised version of the Lame operator (|1U|) . The 
coefficients bk = bk(s) of the corresponding spectral polynomial 

2s 

R2s+i(E) = ]J(E - Ei) = E 2s+1 + hE 2s + b 2 E 2s - 1 + .. + b k E 2s - k+1 + ... + b 2s+1 



up to a sign are the elementary symmetric functions of the eigenvalues: bk — 



(-l) fc e fc , where e x = J2 E h e 2 = J2i<j E i E j> e 3 = T,i<i<k E i E J E k- The ele " 



mentary symmetric functions are related to power sums £?2fc+i(s) = ^2 E i by the 
following well-known relations: 

k 

ke k = y^^-iy~ 1 B 2 j+ie k -j 

3 = 1 

with eo = &o = 1 (see e.g. This implies the following 

Theorem 6. The coefficients b k of the Lame spectral polynomial R 2s+ i(E) are 
related to the elliptic Bernoulli polynomials B 2 j+\{s) by the recurrent relations 

k 



h = -^J2B 2j+1 {s)b k ^. 



3 = 1 

The coefficient b k is a polynomial in s, g±, g 2 , g% with rational coefficients. As a 

2~ 



polynomial in s it has degree 3k and is divisible by (s + l)s(s — l)...(s — [^y^]) 



One can apply this result also to the case of half-integer spin s: in that case all the 
roots of the polynomial R 2s +i(E) are double and correspond to the doubly-periodic 
solutions of the Lame equation. 
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In the reduced case {g\ — 0) the degree of bk drops to [4r] (for k > 1). Using 
the explicit form of the elliptic Bernoulli polynomials given in the Appendix one 
can find the first seven coefficients bk, which in reduced case are: b± = 0, 

62 = ~ uo s{s + 1)(2s ~ 1)(2s + 1)(2s + 3) 

&3 = —Bks(s + l)(2s - 3)(2s - l)(2s + l)(2s + 3)(2s + 5) 
840 

b 4 = 92 s(s - 1) (s + 1) (2s - 1) (2s + 1) (2s + 3) (56s 4 + 76s 3 - 94s 2 + 201s + 630) 
201600 

65 = + 1 ,!f nn s(s-l)(s+l) (2a-3) (2s- 1) (2s+l) (2s+3) (2s+5) (88s 4 +68s 3 -302s 2 +663s+1890) 
1108800 

66 = 201801600 (S ~ 2)(S ~ 1)S(S + 1)(2S ~ 3)(2S ~ 1)(2S + 1)(2S + 3)(2S + 5)X 
(4576s 5 + 12944s 4 - 20720s 3 + 48312s 2 + 597150s + 779625)- 

10378368000 ( * ~ 2)( * ~ 1)s(s + 1)(2s ~ 5)(2s " 3)(2s " 1)(2s + 1)(2s + 3)>< 
(16016s 6 + 89232s 5 + 197160s 4 + 544280s 3 + 2033829s 2 + 3858813s + 2619540) 

2 

, = 9293 (s-3)(s-2)(s-l)s(s+l)(2s-5)(2s-3)(2s-l)(2s+l)(2s+3)(2s+5)x 

24216192000 V A A A A ~ A ~ ) 

(32032s 6 + 189072s 5 + 463440s 4 + 1682920s 3 + 7301418s 2 + 15249213s + 11351340) 

6. Concluding remarks. 

We have shown that for any given k the coefficient bk{s) of the spectral Lame 
polynomial J?2s+i can be computed effectively for all values of parameter s. In 
particular, for fixed s it gives an alternative way to compute the whole polynomial. 
It would be interesting to compare this approach with the classical one going back 
to Halphen and Hermite j2H] further developed recently by Belokolos and Enolski 
P] and Takemura |231 [2*4] following the work of Krichever [J2] . 

However we believe that the elliptic Bernoulli polynomials are of interest by 
themselves. In particular one can expect interesting relations with the arithmetic 
of the corresponding elliptic curves and the representation theory. In this relation 
we would like to mention the elliptic generalisation of the Bernoulli numbers - the 
so-called Bernoulli- Hurwitz numbers BH2k, whose arithmetic was investigated in 

[in eu. 

Another interesting possible relation is with the zeta- function (h(z) = trH~ z of 
the quantum top and its special values. A lemniscatic case 0,3 = ai + a ' 2 could be 
particularly interesting from the arithmetic point of view. 

Recall that the parameter s was originally integer or half- integer (spin). A natu- 
ral question is the role of these values in the theory of elliptic Bernoulli polynomials. 
We conjecture that like in the case of the usual Bernoulli polynomials (see e.g. 26 ) 
these values are the asymptotic positions of the real roots of the polynomials S2/C+1 
for large fc. More precisely, we conjecture that for real s in the bounded interval 
the ratio 

£>2/c+i(s) sin27rs 

as k tends to infinity. Actually, we believe that this is true for each component of 
B2k+i> which is a coefficient at monomial g\g\g\- 
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It is interesting to look at the graphs. In the Fig 1 we show the graphs of the 
coefficients of the polynomial Bis(s) at: (a) g\, (b) gfgi, (c) 313233, and (d) 3I33. 
We normalise each polynomial by dividing it by its first derivative at zero and then 
multiplying it by 2tt. The sinusoidal behavior for small s looks quite plausible. 




Figure 1. 

We would like to mention that the even Bernoulli polynomials (or more precisely 
closely related Faulhaber polynomials) also have elliptic versions related to the 
Lame operator. They were introduced in our recent paper [5j motivated by [7] as 
certain complete elliptic integrals of second kind and have quite different properties. 
The fact that the theory of the Lame equation leads to two different classes of 
polynomials, both related to Bernoulli polynomials (one to odd, another to even) 
seems to be remarkable. To make the picture even more intriguing we would like 
to mention that the integrals in the definition of the elliptic Faulhaber polynomials 
are coming from the formal expansion of the trace of the resolvent of the Lame 
operator (cf. our formula J3Jl). 

Another interesting problem is to investigate the analogues of elliptic Bernoulli 
polynomials related to Sklyanin algebra. It is known after Krichever and Zabrodin 
[H| that Sklyanin's representation [23 gives certain difference analogue of the Lame 
equation, so one can consider the traces of powers of the corresponding generator 
So as functions of the corresponding spin. We would like to mention in this relation 
a very interesting paper .3 , where the traces on the Sklyanin algebra are discussed. 
In particular, the formulas (2.20), (2.21) from this paper give an explicit expression 
of the traces of So an d Sq in terms of elliptic functions, which show that they are 
not polynomials anymore. 

Finally, one can consider our results from the general point of view of the quan- 
tisation of integrable systems. Usually one can find the spectrum in a closed form 
only if the classical system is integrable in elementary functions. The Euler top is 
probably the most natural classical problem integrable in elliptic functions. The 
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question about the nature of its integrability in the quantum case seems to be not 
as easy as it may look. We hope that our paper adds something in this direction 
as well. 
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8. Appendix - The first 8 elliptic Bernoulli polynomials 



B x = 2s +1 



B 3 
B 5 



T2 9i 



+ 



2^0 
60 



s(s + l)(2s + l) 

5 i 2 s(s + l)(2s + l)(3s 2 + 3s-l) 
52 s(s + l)(2*-l)(2* + l)(2* + 3) 



5 i 3 s(s + l)(2s + l)(3s 4 + 6s 3 -3s + l) 
+ Tjk 3152 s(s + l)(2s-l)(2s + l)(2s + 3)(6s 2 + 6s-5) 
+ m 93 s(s + l)(2s-3)(2s-l)(2s + l)(2s + 3)(2s + 5) 



gi 4 s(s + 1)(1 + 2s)(5s 6 + 15s 5 + 5s 4 - 15s 3 - s 2 + 9s -3) 



Pi 52 s(s + l)(2s-l)(2s + l)(2s + 3)(5s 4 + 10s 3 -5s 2 -10s + 7) 
+ «|n 3i33 s 2 (s + l) 2 (2s-3)(2s-l)(2s + l)(2s + 3)(2s + 5) 



+ iSo 32 2 s(s + l)(2s-l)(2s + l)(2s + 3)(4s 4 + 8s 3 -lls 2 -15s + 21) 



33^92 gi 5 s(s + l)(2s + l)(s 2 + s - l)(3s 6 + 9s 5 + 2s 4 - lis 3 + 3s 2 + 10s - 5) 

+ wps 5 l3 5 2 s(s + l)(2s - l)(2s + l)(2s + 3)(20s 6 + 60s 5 -10s 4 - 120s 3 + 44s 2 + 114s -75) 

6 n = + 2^68 3i 2 <73 s(s + l)(2s-3)(2s-l)(2s + l)(2s + 3)(2s + 5)(10s 4 + 20s 3 -4s 2 -14s + 21) 

+ 295^8 S 1 ^ 2 s(s + l)(2s - l)(2s + l)(2s + 3)(40s 6 + 120s 5 -86s 4 - 372s 3 + 242s 2 + 448s -315) 

+ 7m 9 2 9 3 + l)(2s -3)(2s - l)(2s + l)(2s + 3)(2s + 5)(8s 4 + 16s 3 -34s 2 -42s + 63) 
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Bl3 = 



+ 5125120 91*92 

1 3 

2842840 9 1 93 

, 1 2 2 

+ 2562560 ^ 1 92 



s(s + l)(2s + l)(105s 10 + 525s 9 + 525s 8 - 1050s 7 - 1190s 6 
+2310s 5 + 1420s 4 - 3285s 3 - 287s 2 + 2073s - 691) 



s(s + l)(2s - l)(2s + l)(2s + 3) 

(525s 8 + 2100s 7 + 350s 6 - 6300s 5 - 70s 4 + 12810s 3 - 4105s 2 - 11910s + 7601) 
s(s + l)(2s - 3) (2s - l)(2s + l)(2s + 3) (2s + 5) 



(350s 6 + 1050s J 



100s 4 



1950s 3 + 1433s 2 + 2583s - 1650) 



s(s + l)(2s - l)(2s + l)(2s + 3) (1400s 8 + 5600s 7 - 1450s 6 - 23950s 5 
+5438s 4 + 57326s 3 - 24627s 2 - 58215s + 41481) 



ffi<7233 s(s + l)(2s - 3)(2s - l)(2s + l)(2s + 3) (2s + 5) 

(200s6 + 600s5 - 670s 4 - 2340s 3 + 1922s 2 + 3192s - 2475) 



+ 



92 3 s(s + l)(2s - l)(2s + l)(2s + 3) (400s 8 + 1600s 7 - 1640s 6 - 10520s 5 
+8193s 4 + 35786s 3 - 28282s 2 - 48195s + 43659) 

S3 2 s(s + l)(2s - 3) (2s - l)(2s + l)(2s + 3) (2s + 5) 

(80s 6 + 240s 5 - 840s 4 - 2080s 3 + 4401s 2 + 5481s - 7425) 



Bib 



+ 



+ 



+ 



+ 



+ 



5i 



5i 5 92 



1 4 

13178880 9 1 93 



s(s + l)(2s + l)(3s L 2 + 18s 1 ! + 24s 10 - 45s 9 - 81s 8 + 144s 7 + 182s 6 



-345s J 



217s 4 + 498s 3 



- 44s 2 - 315s + 105) 



' 13178880 51 52 



s(s + l)(2s - l)(2s + l)(2s + 3)(42s 10 + 210s 9 + 105s 8 - 840s 7 
-364s 6 + 2730s 5 + 205s 4 - 5540s 3 + 1650s 2 + 5078s - 3185) 

s(s + l)(2s - 3) (2s - l)(2s + l)(2s + 3) (2s + 5)(315s 8 + 1260s 7 
+140s 6 - 3990s 5 + 1265s 4 + 10650s 3 - 5152s 2 - 11352s + 9009) 

s(s + l)(2s - l)(2s + l)(2s + 3(2520s 1 + 12600s 9 + 1750s 8 - 68600s 7 ) 
-13130s 6 + 253630s 5 - 14558s 4 - 557066s 3 + 206601s 2 + 542619s - 360360) 



gi 2 g293 s(s + l)(2s - 3) (2s - l)(2s + l)(2s + 3) (2s + 5)(280s 8 + 1120s 7 - 670s 6 
-5930s 5 + 3047s 4 + 17284s 3 - 11237s 2 - 21054s + 18018) 



5i52" 



5i53 



5253 



s(s + l)(2s - l)(2s + l)(2s + 3)(1120s 10 + 5600s 9 - 2400s 8 - 43200s 7 
-8814s 6 + 201162s 5 - 60127s 4 - 517124s 3 + 256797s 2 + 557766s - 405405) 

s 2 (s + l) 2 (2s - 3)(2s - l)(2s + l)(2s + 3) (2s + 5) 

(80s 6 + 240s 5 - 840s 4 - 2080s 3 + 4401s 2 + 5481s - 7425) 

s(s + l)(2s - 3) (2s - l)(2s + l)(2s + 3) (2s + 5)(80s 8 + 320s 7 
-600s 6 - 2920s 5 + 4037s 4 + 13314s 3 - 16959s 2 - 24156s + 27027) 
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